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A phenomenological diffusion law L(t) oc i' 3 , where L(t) measures the spreading of a wave-packet in 
a time t, is assumed for perfect quasicrystals. We show that it affects their conductivity with striking 
differences compared to the case of periodic metals. In the absence of defects the dissipative part 
of conductivity is non zero even at low frequencies contrary to the case of crystals. Also if j3 < 1/2 
the d.c. conductivity increases when disorder increases and the so-called Drude peak, characteristic 
of metals, is replaced by a dip . Experimental results are briefly discussed. 



PACS numbers: 71.23. Ft, 78.20.-e , 72.15.-v 

Since the pioneering work of Kohmoto et al [3] nu- 
merous studies of quasiperiodic Hamiltonians have given 
evidences that the associated eigenstates are multifractal 
and critical with an algebraic decrease of their envelop 
at large distances || ||]. These states lead to a peculiar 
quantum diffusion which is neither ballistic as in crys- 
tals nor diffusive as in metallic disordered systems but 
follows at sufficiently large time a power law L(t) oc t@ 
where L(t) measures the spreading of a wave-packet in a 
time t . (3 depends on the energy of the wave-packet 
and on the Hamiltonian but is restricted to < j3 < 1. 
For realistic models, periodic approximants of quasicrys- 
tals |-0 have narrow bands, in agreement with the 
tendency to localization. For some phases these com- 
putations have even confirmed a multifractal character 
of eigenstates Q . Electronic transport has been studied 
through Bloch-Boltzmann type approximation assuming 
a peculiar bandstructure ]9|-pd| and this provides inter- 
esting information. Yet this approach does not consider 
the anomalous diffusion. 

This paper presents a formalism that allows to discuss, 
for the first time to our knowledge, how the a.c. con- 
ductivity is affected by the anomalous diffusion. For d.c. 
conductivity of quasiperiodic systems with some disor- 
der, the role of anomalous diffusion has been considered 
by several authors []l^-|l4j|. Using the Einstein formula 
a = 2e 2 n(Ep)D(Ep) where n(Ep) is the density of states 
at the Fermi energy per spin and D{Ep) the diffusivity 
of states at the Fermi energy one expects that the zero 
frequency conductivity is : 



cr DC (r) ~ 2e 2 n{E)Ar 



2/3-1 



(f) 



where r is the time beyond which the propagation be- 
comes diffusive due to disorder. In this scheme the prop- 
agation on time scale t > t is diffusive with D(Ep) — 

L ^ — At 2 ^' 1 . For (3 < 1/2 the conductivity increases 
with disorder as observed in some quasicrystalline alloys 
p5[ . From a dimensional argument one expects that the 
frequency dependent conductivity dAci 1 ^) = £" \dc( t = 
^j) at the lowest frequencies where B is some 

complex number with a modulus of order unity . The 
real part of B is related to energy dissipation by the sys- 



tem. For crystals it is zero but the dimensional argument 
alone cannot give its value in quasicrystals. 
In this paper we study the zero or low frequency con- 
ductivity of a quasiperiodic Hamiltonian of independent 
electrons with and without static disorder. It is shown 
first that a suitably defined average of the conductiv- 
ity at frequency to is related to a quantity cr(w) which 
has the analytical properties of a conductivity and is 
directly related to quantum diffusion. The behaviour 
of <j{uj) with frequency or disorder is thus representa- 
tive of that of the conductivity and is studied in de- 
tail. Assuming an anomalous diffusion law in the per- 
fect quasicrystal we obtain an exact result for the con- 
ductivity <7p ure (w) of the quasiperiodic system without 
defects. Then we consider the role of static disorder and 
introduce a relaxation time approximation for the ve- 
locity auto-correlation function which physical meaning 
and limits are discussed. This approximation allows to 
derive a generalized Drude formula for ct(cj) that inter- 
polates between the above <jdc{ t ) (JlJ) and crp ure (uj) as 
a function of the parameter lot. This formula implies 
a drastic deviation of optical conductivity from the case 
of crystals. We comment shortly on the conductivity of 
several quasicrystalline phases. 

Let us define several quantities and give some exact rela- 
tions between these quantities. We define first < f >e 7 
the average diagonal element of an operator / among 
states at energy E : 



< f >E-' 



Tr[6(E - H)f] 
Tr[5{E - H)] 



(2) 



Tr means a trace of an operator and 8(E — H) projects 
on eigenstates of the Hamiltonian H at energy E. Next 
we define a measure X 2 (E,t) of the spreading of wave- 
packets at energy E along the x-direction : 



X 2 (E,t) =< (X(t)-X(0)) 2 >! 



(3) 



where X{t) is the position operator along the x- 
direction at time t in the Heisenberg representation. 
X 2 (E,t) is a real and pair function of time. We define 



1 



also the auto-correlation function C(E, t — t') for the ve- 
locity V x (t) — ^Jf* between t and t' and which from (|^) 
satisfies: 



C(E,t- =< V x (t)V x {t') + V x (t')V x (t) >! 



dX 2 (E, t) 
dt 



dt'C{E,t') 



(4) 



(5) 



C(E, t) is a real and pair function of t. For pure qua- 
sicrystals the phcnomcnological diffusion law reads : 



Xl(E,t)~At 2p 



(6) 



Co(E, t) = d2X ^ ,t] * A2[3(2[3 - l)t 2 ^ 2 (7) 

(0) is obtained from (0) and (^). In periodic approxi- 
mants, the scaling law Xq(E, t) ~ Ai 2/3 will be valid only 
up to some maximum time ti beyond which the period- 
icty is effective. Then quantities deduced from the scal- 
ing law (see below) will be valid for time scales r < t<z or 
frequencies scales co > j-. The Guarneri inequality fl6| l 
implies a singular continuous density of states of a three 
dimensional quasicrystal for /3 < 1/3. Yet if we consider 
periodic approximants we can assume that the density 
of states is absolutely continuous even for f3 < 1/3 but 
with the above restriction on time and frequency scales. 
From its definition (|) C(E,t = 0) = 2 < V 2 > E is fi- 
nite which implies that the propagation is ballistic at the 
lowest times. Thus the anomalous diffusion law can be 
valid only at sufficiently large times t > t± . It occurs only 
when the wave-packet has sufficient time to " feel" specific 
atomic arrangements |T7j , |l8fl and longer range quasiperi- 
odic order. This remark is important for the sum-rule 
satisfied by conductivity (|l5|). 

We introduce F{E,uS) that appears in the formula JToj ] 
for the component a XjX (fj,,u) of the conductivity tensor 
and satisfies: 

F(E, u) = ——Tr[5(E - H)V X S(E + huj- H)V X ] (8) 



S{E, w) = F(E, lo) + F(E - hu, w) 



F(E,w) in terms of the velocity correlation function. 
This relation is often valid in disordered systems at low 
frequency j2lj|. Yet in quasicrystals there are important 
variations of electronic structure on small energy scale 
(tenths or hundredth of electron Volt and possibly less 
P,P|) and the above assumption could have a vanishingly 
small range of validity . Since we are interested in gen- 
eral trends we consider an average of the conductivity at 
frequency u over a range of chemical potential values (we 
could also consider the conductivity at very high temper- 
ature kT >> hu). More precisely we define averages of 
functions g(E) by a convolution with a positive function 
P(x) such that / dxP(x) = 1 : 



< g{E) > P = J dxP{x)g{E - x) 



(11) 



One can choose for example P(x) — 2W ex P(~w^ 
where W is its width. From the Kubo-Greenwood for- 
mula (H at T = : 

A* 

/dE 
— <F(E,uj)>p (12) 

where Re means the real part. If W » hu F(E,u) 
and F(E — hu,u) have nearly the same weight in < 
F(E,u) > P dH). Then only their sum S(E,w) is rel- 
evant and F(E,u) can be replaced by s ^' u ^ in (p^). 
One can write : 

f dE 

< Rea XlX ((i,u) > P ~ / — < Red(E,u) > P (13) 



fi—hu) 



cr(/i, uj) = e n(fj>) / exp{iwt)C(iJLit)dt 



(9) 



(14) 



since for u real S(E,u) = 2Rea(E,u).The relative 
difference between the two members of ( |13| ) is bounded 
by a term of order ^ for u real [p0|] . a(fj,,ui) has the 
dimension of a conductivity. It is an analytical function 
of u in the upper half of the complex plane. For large 
u ct(/z, u) oc — and the Kramers-Kronig relations can be 
proved. Finally the usual sum rule is also valid : 



S(E,u) = e 2 n(E) / exp(iut)C(E,t)dt 



(10) 



51 is the volume of the system. Starting from (g), 
( |To| ) is easily shown by using the Fourier expansion of 
S(E — H) and the cyclic property of a trace pC[ |. If 
F(E, u) depends wealky on the energy E at a given u 
then F(E,u) ~ F(E - hu,u) ~ S(E,u)/2 which gives 



Rea(fj,,u)du = 7re 2 n(^) 



C(n,t = 0) 



7re 



2m* 



15) 



where the last equality holds for free electron like systems 
with electron density n and effective mass m*. From now 
on we focus on the properties of a(fi, u) since ( |l3| ) shows 
that its variation with frequency or with defects drives 
the variation of the average conductivity < a(p,u) >p. 
Due to Kramers-Kronig relations we expect also that the 
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imaginary part of a(/j,, ui) is representative of that of the 
averaged conductivity < a(fj,,u>) >p as discussed in [ p0[ . 
Using the velocity correlation ([?]) associated to the 
anomalous diffusion law (|^) we deduce ap ure (^,,Lo) : 



a Pure (n,oj) ~ e z n{ii) AT {2(3 +!)(-) 



2(3-1 



(16) 



since the anomalous diffusion is valid only at t > t\ the 
above formula is valid for w < j-. T(z) is the Gamma 
function and is of order one for < /3 < 1 |2^] . 
The real part of this expression is strictly positive if 
< (5 < 1 since i?e(i 2/3 - 1 ) = cos(§(2/3 - 1)) . Thus 
there is dissipation even for a perfect quasiperiodic sys- 
tem. This is in contrast with the case of crystals ((3=1) 
for which the dissipation is zero. A natural interpretation 
is obtained by comparing the bandstructure of crystals 
and that of approximants of quasicrystals in the limit of 
infinite unit cell. In a crystal the number of bands is 
finite and when lo tends to zero there is no more possibil- 
ity for vertical interband transition at energy Hlo between 
an occupied state and an empty one. Since the system 
cannot absorb energy the real part of the conductivity is 
zero. On the contrary the number of bands of approxi- 
mants with large unit cells increases indefinitely with the 
size of the unit cell. In the limit of infinitely large unit 
cell interband transition are possible even at the small- 
est frequencies and the dissipative part of conductivity is 
non zero. 

We consider now the effect of static disorder. For weak 
static disorder numerical studies [^3| give evidences that 
the propagation is diffusive on time scale t > r. There 
could be exceptions [^4f at very strong quasiperiodic po- 
tential for which the quantum diffusion stays anomalous 
even with static disorder. In this case the result ( |l6| ) is 
still valid. Our aim is to find an interpolation formula 
from the low frequency regime lot << 1 (|j) to the high 
frequency regime lot >> 1 ([l6]) . Let us assume a relax- 
ation time approximation for C(E,t) : 



C(E,t)~C {E,t)exp(-^) 

T 



(17) 



where Co (E, t) is the velocity correlation function of 
the system without disorder. The physical meaning is 
that due to disorder the electron is scattered in states 
with no or smaller velocity correlation leading to a damp- 
ing of the velocity correlation function. Using (|l7]) we see 

that in (|^) dx ^ is constant, and the propagation is 
diffusive, for t > r. The approximation ( |l7| ) could be 
generalized by replacing the exponential by an average 
of damped exponentials. It will be clear that due to the 
linearity of equations this does not change our qualitative 
conclusions. 

Assuming ©, using X 2 (E, 0) = (see®), = 

at t = (see (||)) and defining z = {- — ico), <r(/i, lo) Jl4| ) 
is equivalent to : 



a(n,co) ~ e 2 n(fj,)Az 2 / dtXQ(E,t)exp(—zt 



d{Li,Lo)~e 2 n{n)AY{2l3 + l){- — ) 

1 — VjJT 



2/3-1 



(18) 



(19) 



( |18| ) expresses a direct relation between quantum diffu- 
sion and optical conductivity. (|l£) is obtained from (IS) 
and (^|). Since the diffusion law (^3) is valid for t > t\ (1£) 
applies only for lo < ~ and r > t\ and this generalized 
Drude formula does not satisfiy the sum rule (|l5|). The 
point is that it interpolates between ( |l6| ) for lot » 1 
and (Q) for lot << 1 (up to a numerical factor of order 
one). 
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FIG. 1. Variation of R(lo) = 
for several values of a = (2/3 



i?ecr(/i, lo 
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as a function of lut 



Figure (|i|) shows R{lo) = ith I as a function of 

cr(/x,0) 

lot. For a = 2/3 — 1 > the so-called Drude peak is 
less marked than in crystals yet the distinction may be 
difficult to establish experimentally. The deviation from 
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the case of crystals is drastic if a = 2/3 — 1 < : instead 
of a Drude peak there is a plateau or a dip. 
It is important to note that due to quantum interfer- 
ences similar to those of disordered systems jL2| the veloc- 
ity correlation function should become negative on time 
scale t > t. This effect is beyond the relaxation time 
approximation (JlJ) and becomes quantitatively impor- 
tant close to the metal-insulator transition |2l| . Accord- 
ing to the scaling theory of localization the parameter 
g = <JDrudeXo(T) must be compared to the universal 
critical conductance for three dimensional systems g c . If 
g >> g c the corrections are negligible and the relaxation 
time approximation is essentially correct. If g is compa- 
rable to g c the relaxation time approximation is incorrect 
and if g < g c the quantum interferences even lead to an 
insulating state. Quantum interferences are effective on 
time scale t > t or oj < i and it has been argued 0] 
that the system behave like a standard disordered sys- 
tem at these time scales. In any case one expects that 
( |l9| ) is correct for j- > u> > \. Even though the formal- 
ism is for static disorder, we believe that the generalized 
Drude formula (19) is better for an inelastic scattering 
since then quantum interferences cannot appear. 
For many decagonal phases d-AlCuCo and icosahedral 
phases i-AlCuFe it seems that ~ 310 _2 eV for I 13 

T < 400K H) then the full range in Figure ([!]) is 
hcu ~ 0.6eV\ A detailed comparison of the present theory 
with experiment is difficult since we consider only aver- 
ages of the conductivity and also because quantum inter- 
ferences effects EB] that are not included in the relaxation 
time approximation should play a role at low frequency 
(W(T) < 1). Yet since i-AlCuFe §| and i-AlPdRe @ 
do not exhibit a Drude peak they are similar to the case 
2/3 — 1 negative. Correlatively (in agreement with (|ig|)) 
the d.c. conductivity increases with temperature and 
probably also with the addition of defects. For d-AlCuCo 
p7| , which consists of quasiperiodic planes stacked peri- 
odically in the z-direction, the dissipative part of conduc- 
tivity along the quasiperiodic planes depends weakly on 
frequency which is consistent with 2/3—1 close to zero. 
Correlatively the conductivity in the quasiperiodic planes 
depends weakly on temperature. Along the z direction 
the conductivity of d-AlCuCo is that of good metals. 
To conclude let us recall that a suitably defined average of 
the conductivity can be related to the quantum diffusion 
in the system X (fi,t) as expressed by (O, |l8|). Strik- 
ing differences with crystals occur when Xo(fi, t) deviates 
from a ballistic law and this seems to be in agreement 
with experimental results. 
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